Massive stars in the CI 1813-178 Cluster. An episode of massive star formation 

in the W33 complex. 



Maria Messineo 1,2 , Ben Davies 3,5 ' 10 , Donald F. Figer 3 , R.P. Kudritzki 7 ' 9 , Elena Valenti 4 , 
Christine Trombley 3 , F. Najarro 8 , R. Michael Rich 6 

messineoOmpif r-bonn . mpg . de 
ABSTRACT 

Young massive (M > 10 4 M Q ) stellar clusters are a good laboratory to study the 
evolution of massive stars. Only a dozen of such clusters are known in the Galaxy. Here 
we report about a new young massive stellar cluster in the Milky Way. Near-infrared 
medium-resolution spectroscopy with UIST on the UKIRT telescope and NIRSPEC on 
the Keck telescope, and X-ray observations with the Chandra and XMM satellites, of 
the CI 1813-178 cluster confirm a large number of massive stars. We detected 1 red 
supergiant, 2 Wolf-Rayet stars, 1 candidate luminous blue variable, 2 Olf, and 19 OB 
stars. Among the latter, twelve are likely supergiants, four giants, and the faintest three 
dwarf stars. We detected post-main sequence stars with masses between 25 and 100 
Mq. A population with age of 4-4.5 Myr and a mass of ~ 10000 M can reproduce such 
a mixture of massive evolved stars. This massive stellar cluster is the first detection 
of a cluster in the W33 complex. Six supernova remnants and several other candidate 
clusters are found in the direction of the same complex. 

Subject headings: stars: evolution — infrared: stars 
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Introduction 



An understanding of the mechanisms of formation, evolution, and end state of massive stars 
is fundamental for the studies of galaxies at all redshifts. Massive stars contribute to the chemical 
enrichment of the interstellar medium with their strong winds and by exploding as supernovas. 
Massive stars are the most luminous stars, can easily be detected in external galaxies, and provide 
distance estima tes. They are the source s of the most energetic phenomena in the Universe, gamma 
ray bursts (e.g. IWooslev Bloom! 12000 ). 



The availability of large surveys of the Galactic plane at radio and infrared wavelengths opens a 
golden epoch for studying the formation, evolution, and environments of massive stars. More than 
1500 new candidate stellar clusters h ave been discovered, and among them several young clusters 
rich in massive stars may be hidden (jMessineo et al.l l2009). 



In iMessineo et al.l (|2008l ) (hereafter referred to as Paperl), we presented the serendipitous 
discovery of a young massive cluster, CI 1813-178, in the Galactic disk at 1=12°, with a spectro- 
scopically identified population of massive stars, a red supergiant star (RSG), two blue supergiants 
(BSG), and one Wolf-Rayet (WR) star. Here, we present a follow-up study of the cluster. Near- 
infrared photometry and spectroscopy, and X-ray observations, reveal a large number of massive 
stars. 

The cluster is located in the W33 complex, and is associated with two supernova remnants 
(SNR), SNR G12. 72-0.00 and G12. 82-0.02, and the highly magnetized pulsar associated with the 
TeV 7-ray source HESS J1813— 178. Interestingly, the W33 complex appears to contain several 
other candidate ste llar clusters, and several SNRs. Clusters do form in large complexes (e.g. 
Beuther et al.l 120071 ). and their spatial distribution varies from cloud to clo ud, indicating that 
several external and internal triggers can be at work (e.g. IClark et al.l 120091 ) . W33 is an ideal 
laboratory to investigate various issues about massive stars and multi-seeded star formation, and 
to verify the presence of triggered sequential star formation, which is suggested by the presence of 
SNRs. The association of the stellar cluster with two SNRs can shed light on the initial masses of 
the supernova progenitors, and on the fate of massive stars. 

We describe the observations and data reduction in Sect. 2, and we analyze the spectra and 
the cluster color-magnitude diagram (CMD) in Sect. 3. A discussion on the massive members and 
spectrophotometric distances is presented in Sects. 4 and 5. Cluster age and mass are derived in 
Sects. 6 and 7. An overview of candidate clusters in the direction of the W33 complex is given in 
Sect. 8. Finally, our findings are summarized in Sect. 9. 
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Observations and data reduction 



2.1. IR photometry 

Pho tometric measuremen ts from the near-infrared Two Micron All Sky Survey (2MASS) 
catalog (jSkrutskie et al.l 120061) . the Galact ic Legacy Infrared Mid-Plane Survey Extraordinaire 
(GLIMPSE) catalog ([Benjamin et al.ll2003r). as wel l as from the Naval Observatory Merged As- 
trometric Dataset (NOMAD) (jZacharias et al.l 12004 ) were used. 



Images of the central cluster region with the J, H and K s filters were obtained during two 
nights of observation, 2 008 June 23-240, using the near-IR camera SofI mounted on the ESO NTT 
(|Moorwood et al.lll998l ). We used SofI in large-field mode with a pixel size of 0.288", and a total 
field of view of 4.9' x 4.9'. For each filter, we performed a random dithering pattern of 14 images, 
and reached a total exposure time of 840s, 1176s, and 1680s in J, H and K s , respectively. Each 
image is a combination of 50, 70, and 100 exposures, each one 1.2 s long, in J, H, and K s bands, 
respectively. 

Data reduction was performed using standard IRAF routine. For each filter, we obtained a sky 
image by median-combining the dithered frames, and we subtracted the sky image from each frame, 
Flat fielding was performed by using the SpecialDomeFlat template, which applies the appropriate 
illumination correction, as described in the SofI User Manual. Finally, the dithered frames were 
averaged into a single image. Standard crowded-field photometry, including point-spread function 
modeling, was carried out on each image using DAOPHOT II/ALLSTAR (Stetson 1987). The 
internal photometric accuracy was estimated from the rms frame-to-frame scatter of multiple star 
measurements (0.03 mag< o j^(Jh^ok 3 <0.06 mag for 8< J~H~K S <18). The instrumental 
magnitudes were converted into the 2MASS photometric systeiro The 2MASS catalog was used 
as an aa.rome.ric reference framfi For stars that were saturated on the SofI images, magnitudes 
were taken from the 2MASS catalog. 



2.2. Near-IR spectroscopy 

Spectroscopic observations of the brightest stars were carried out with NIRSPEC at the Keck 
Observatorjfl under program H243NS (PI: Kudritzki) on 2008 July 24. The if-filter and a 42" x 
0.570" slit were used, covering from 2.02 to 2.45 //m with a resolution of R=1700. For each target, 

1 Based on data taken within the observing program 081.D-0371. 

2 An overall uncertainty of ±0.05 mag in the zero-point calibration in all the three bands has been estimated 

3 The astrometric procedure provided rms residuals of ~0.2" in both the right ascension and declination. 

4 Data presented herein were obtained at the W.M. Keck Observatory, which is operated as a scientific partnership 
among the California Institute of Technology, the University of California and the National Aeronautics and Space 
Administration. The Observatory was made possible by the generous financial support of the W.M. Keck Foundation. 
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Fig. 1. — 2MASS i^-band image of the cluster. The big box shows the location of the field observed 
with Sofl. Stars that were spectroscopically observed and are likely cluster members are marked 
with small boxes. Non-member stars are marked with triangles. 
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two nodded exposures of 10s each were taken. We observed a total of 23 stars with NIRSPEC (see 
Tables ffland ED. 

Additional spectroscopic data were taken with the UKIRT 1-5 micron Imager Spectrometer 
(UIST) at the UKIRT Observatory |. We used the short-K grism in combination with a 120" x 
0.12" slit, covering from 2.00 to 2.26 fim at a resolution R=1800. We also used the long-K grism 
to cover the spectral range from 2.204 to 2.513 /im at a resolution of R=1900. Typically, each 
target was observed with the long-K grism. When CO bands were not detected, a second spectrum 
was taken with the short-K grism. Integration times varied from 30 to 60 s per exposure, and the 
number of exposures varied from 8 to 20 s. The 39 stars observed with UKIRT, including a few 
chance detections, are listed in Tables Q] and [71 

Pairs of frames with nodded positions were subtracted and flat-fielded. The stellar traces were 
straightened using a two dimensional de-warping procedure. We wavelength calibrated the spectra 
with arc lines. We corrected the spectra for atmospheric absorption and instrumental response 
by dividing the observed spectra by the spectrum of a reference star. Reference stars were of B 
type (from B2 to B9). The Br7 and Hel lines of the telluric spectrum were eliminated with linear 
interpolation. 

The signal to noise of the spectra varied from 40 to 150. 



2.3. X-ray data 



The cluster is located in the vicinity of the HESS J1813-178 pulsar wind nebula (jHelfand et al 



20071 ) . Several X-ray observations o f this region ha v e bee n performed in order to identify the pulsar 
associated with the HESS source. iHelfand et al.l (|2007l ) presented a catalog of 75 point sources 
detected with the Chandra satellite. We cross-correlated the Chandra catalog with the 2MASS 
point source cata log, and identified a total of 44 matches (Paperl). Using the XMM satellite, 
Funk et al.1 (120071 ) detected seven X-ray sources in the surrounding area of HESSJ1813-178. 



3. Analysis 

3.1. Spectral classification 

Spectral class i fication was performed by comparing the spectra with spectral atlases (e.g. 

Hanson et allll996l . bood : iFiger et al.lll997i : iMartins et aDl2007i : kleinmann & Halllll986l : IWallace & Hinkk 
1996h. 



5 The United Kingdom Infrared Telescope (UKIRT) is operated by the Joint Astronomy Centre on behalf of the 
Particle Physics and Astronomy Research Council. 
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We observed a total of 60 stars, and detected 24 early-type stars and 36 late-type stars. In 
addition, we considered the previous observations reported in Paperl, which added an extra early- 
type star (#11) to our new sample. The observed stars were divided into candidate cluster members 
and non members, and are listed in Tables [H and [71 We first listed the sample of stars given in 
the Table 1 of Paperl, using the same identification numbers, then we appended the new targets 
ordered in K s magnitude. The spectra of the newly detected early-type stars are shown in Figs. 
and [31 

In the following, we describe the spectra by grouping the stars into late-types and early-types. 
Early types divided in four subgroups (Wolf-Rayet stars, candidate luminous blue variables, O-type, 
and B-type stars). 



3.1.1. Red Supergiants 

In the -ff-band, the spectra of late-type stars have CO band-heads at 2.29 fim. By measuring 
the equivalent widths of the CO bands, it is possible to distinguish between giants and supergiants 
(see the Appendix). 

Among the detected late-type stars, four stars, #1, #32, #38, and #39, are found to have 
equivalent widths of the CO bands, EW(CO), typical of RSGs. 

Star #1 is a K2-K5I cluster member, as discussed in Paperl on the basis of high-resolution 
data. 

Star #32, #38, and #39 have EW(CO)s typical of RSGs with spectral types M2.5, M3.5 
and Ml, respectively. Their photometric properties indicate that they are unrelated to the stellar 
cluster. We further discuss these stars in the following sections. 



3.1.2. Wolf-Rayet stars 



Star #4 is a known WR star of type WN7 (number #8 in lHadfield et al.ll2007l ). Star #7 has 
similar emission lines at 2.114 /im, 2.166 //m, and 2.189 /im. By comparing it with the spectral 
atlas of lFiger et al.1 (|1997l ). star #7 appears to be another WN7 star (see Fig. (H). 



3.1.3. Candidate Luminous Blue Variables 



The spectrum of star #15 is dominated by He I lines at 2.058 fj,m and 2.112 /im, and Br7 line 
in emission. Faint Mgll lines are detected at 2.138 /im and 2.144 /im. The type of detected lines 
and the line profiles indicate that star #15 is a P-Cygni B supergiant. Its spe ctrum is similar t o 
that of the candidate luminous blue variable (cLBV) star in the W51 region (jClark et al.ll2009l ). 
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More information on the stellar photometric variability is needed to firmly confirm that it is another 
LBV star. 



3.1.4- O-type stars 

Spectra of O and B stars have hydrogen lines, helium lines, and other atomic lines (e.g. CIV, 
NIII, CIII, Mgll). Hydrogen and helium lines are usually seen in absorption; they may be in emission 
in supergiants. The Br7 line is the only hydrogen line detectable in K-band. The Hel transitions 
at 2.058/mi and 2.112 fim are usually detected in late-0 and early-B stars. The Hell line at 2.189 
/jm appears only in O-type stars. With our resolution, the Hell line at 2.189/iin may be present 
down to 09.51 or 09V. 

The spectrum of star #6 has a strong Hel line at 2.058 fim in absorption, CIV lines at 2.069 
/im and 2.079 fim in emission, strong emission line at 2.11 //m (probably due to Hel/NIII/CIII), 
and a Hell absorption line at 2.189 /um. NIII at 2.25/im is seen in emission. There is an indication 
for a faint CIV line at 2.083 fim, and a NIII line at 2.10 /xm. The simultaneous presence of Hel, 
Hell, and CIV lines is typical of an O-type star. The strengths of the Hel and Hell absorption 
lines, as well as the detection of NIII at 2.25/mi, indicate a 06-07f+. 

The spectrum of star #16 has Br7 line in emission, a Hel line at 2.058 fim with P-Cygni 
profile, and a Hel/NIII line at 2.112 /im. A Hell line at 2.189 /im is seen in absorption. A faint 
line is seen in emission at 2.105 /ini. Only an indication for a CI V line at 2.08 um i s visible. The 



spectrum resembles that of a 08-09If star, similar to HD151804 (IMorris et al.lll996l ). 



Stars #18, #19, and #20 have Hel lines at 2.058 and 2.112 //m in absorption, a NIII line 
at 2.115 //m in emission, and Br7 line in absorption. A Hell line at 2.189 [im clearly appears 
in absorption. The absence of CIV lines, and the presence of NIII and Hell lines, suggest a late 
O-type (07-09). 

The spectra of stars #5, #8, and #9 have Hel lines at 2.058 fim, and 2.112 /im in absorption, 
and Br7 line in absorption. There is an indication for a NIII line at 2.115 in emission, and a Hell 
line at 2.189 /im in absorption. Star #5 has a clear detection of a NIII line. The absence of CIV 
lines, and the presence (or indication) of NIII and Hell lines, suggest that this is a late O-type 
(07-09). 



3.1.5. B-type stars 

Here, we list all detected early-type stars, which do not have Hell lines. These may be B-stars 
or 09 dwarfs. With our resolution, the Hell line at 2.189/im may be present down to 09.51 or 09V. 
Once the Hell line disappears, assigning a spectral type and luminosity class can turn into a quite 
degenerate issue. This degeneracy may be partly broken if the absorption components of the Br7 
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line and of the Hel lines are sufficiently filled in by the stellar wind contribution. Typically, a Hel 
line at 2.112//m in absorption is observed down to B8Ia stars, while the same line is not detected 
in dwarfs/giants later than B3V. Such matters are highly dependent on the quality of the spectra. 



From a comparison with the atlas of lHanson et al.1 (119961 ). stars #2, #3, #11, #12, #13 



#14, #17, #21, #22, #23, #24, #25, and #26 appear to be likely stars with types from BO to 
B2. 

The spectra of stars #2 and #13 have a Hel line at 2.058 fim in emission, a He l line at 2.112 
//m, and Br7 line in absorption. From the relative strengths of the Br7 and Hel lines (jHanson et al 



19961 ). the presence of Hel in emission at 2.058/xm, and the lack of Hell lines, these stars appear to 



be B0-B3 supergiants. 

The spectra of stars #3, #22, #24 , #25, and #26 have a Hel line at 2.112 /j>m in absorption, 
as well as Br7 line in absorption. Star #17 has Br7 line in absorption, and a trace of Hel at 
2.112/mi. These are likely B0-B3 stars or 09-B3 stars (depending on the luminosity class). 

The spectra of stars #11, #21, and #23 have Br7 line in absorption, Hel lines at 2.058 and 
2.112 /im in absorption. These spectra are typical of early B-type stars. 

Star #12 has Br7 line in absorption. The lack of Hel al 2.112 in absorption suggests that star 
#12 has a spectral type later than a B8I or B3-B4V star. 

The spectrum of star #14 has Br7 in absorption, a Hel line at 2.112 /jm in absorption, an 
indication for a Hel line in emission at 2.058 /im, and a faint Hel line at 2.113 /jm in emission. The 
lack of Hell at 2.189/im and presence of Hel lines suggest that star #14 is a BSG star with spectral 
type between BO and B3. 

The spectrum of star #11 is shown in Paperl. Stars #2 and #3 were also observed with 
NIRSPEC in April 2008 (Paperl), and we obtained the same spectral classification. 



3.2. Photometric results 

3.2.1. Cluster size 

We assume as a cluster center the flux weighted centroid of a 2MASS if s -band smoothed 
image (RA=18:13:24.15, DEC=-17:53:29.64). We take as a cluster radius the average distance 
from the cluster center where the surface brightness becomes equal to the average brightness of 
the surrounding field (3.5'). The 61 spectroscopically observed stars are located within 6.5' from 
the cluster center, but predominantly (50) within the cluster radius. Within 3.5' from the cluster 
center, we observed 41 stars out of the 48 stars with K s < 9.5 mag, and 30 stars out of the 31 stars 
with K s < 9.0 mag. 

A fraction of 48% of stars brighter than i^ s =9.5 mag, within the cluster radius, have early types, 



Table 1. List of cluster members spectroscopically observed. 



ID 




Ra 


Dec 


B 


V 


R 


J 


H 


K s 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


tel. 


Sp 


l a 


18 


13 


22.26 


- 1 7 54 15 55 


14.74 


13.68 


12.73 


5.46 


4.23 


3.79 


5.07 


4.05 


3.57 




keck 


RSG 


2 


18 


13 


23.26 


-17 53 26 54 


20.08 


15.03 


14.50 


8.64 


7.72 


7.22 


6.94 


6.78 


6.63 


6.68 


keck 


B0-B3 


3 


18 


13 


24.43 


-1 7 52 ^6 75 


20.42 


15.74 


14.68 


9.20 


8.25 


7.79 


7.47 


7.38 


7.30 


7.30 


keck 


B0-B3 


4 


18 


13 


14.19 


-17 53 43 60 


20.48 


17.57 


15.20 


9.62 


8.60 


7.94 


7.34 


6.93 


6.70 


6.35 


keck 


WN7 


5 


18 


13 


23.71 


-17 50 40.41 


19.51 


16.41 


14.94 


9.96 


9.06 


8.56 


8.23 


8.10 


7.99 


8.06 


keck 


07-09 


6 


18 


13 


33.82 


-17 51 48 89 






16.12 


10.29 


9.15 


8.57 


8.00 


7.90 


7.75 


7.75 




06-07f 


7 


18 


13 


22.52 


-17 53 50 14 


18.63 


15.60 


16.57 


10.30 


9.27 


8.66 


8.03 


7.70 


7.49 


7.30 


ukirt 


WN7 


8 


18 


13 


14.49 


-17 54 01 77 


19.91 


15.47 


15.71 


9.86 


9.12 


8.75 


8.43 


8.33 


8.31 


8.31 


ukirt 


07-09 


9 


18 


13 


21.91 


-17 51 35 91 


21.17 




16.24 


10.74 


9.80 


9.34 


8.95 


8.90 


8.83 


8.91 


ukirt 


07-09 


11 


18 


13 


23.62 


-17 53 24 46 








11.74 


10.85 


9.59 


10.16 


10.32 


10.02 


99.99 


keck; 


B0-B3 


12 


18 


13 


41.33 


-17 55 41.14 




17.42 


18.09 


9.07 


7.50 


6.75 


6.41 


6.31 


5.87 


5.95 


keck 


B9-A2 


13 


18 


13 


29.18 


- 1 7 55 14 55 


21.09 




15.66 


8.79 


7.56 


6.85 


6.80 


6.15 


6.03 


5.99 


keck 


B0-B3 


14 


18 


13 


26.47 


-17 52 49.48 


16.27 


14.25 


13.01 


8.21 


7.38 


6.91 


6.81 


6.44 


6.38 


6.32 


keck 


B0-B3 


15 


18 


13 


20.98 


-17 49 46.88 


19.67 




15.40 


8.93 


7.78 


6.96 


6.76 


5.91 


5.71 


5.42 


keck 


cLBV 


16 


18 


13 


19.13 


-17 52 57.43 






16.90 


9.43 


8.10 


7.33 


6.94 


6.47 


6.40 


6.19 


keck 


08-09K 


17 


18 


13 


27.83 


-17 52 33.19 




13.46 


13.13 


9.70 


8.97 


8.52 


8.29 


8.23 


8.16 


8.20 


keck 


B0-B3 


18 


18 


13 


25.25 


-17 53 55.05 


20.13 


16.60 


14.02 


10.05 


9.05 


8.61 


8.25 


8.09 


7.96 


8.04 


ukirt 


07-09 


19 


18 


13 


25.20 


-17 53 55.83 


19.49 


16.28 


16.20 


10.25 


9.25 


8.72 


8.41 


8.25 


8.21 


8.20 


ukirt 


07-09 


20 


18 


13 


21.54 


-17 54 35.65 








10.60 


9.61 


9.14 


8.70 


8.62 


8.49 


8.59 


ukirt 


07-08 


21 


18 


13 


25.06 


-17 53 09.34 




16.84 


13.98 


10.73 


9.78 


9.29 


9.03 


8.86 


8.85 


8.77 


ukirt 


09-B3 


22 


18 


13 


19.25 


-17 53 22.62 








11.42 


10.10 


9.34 


8.76 


8.57 


8.40 


8.50 


ukirt 


09-B3 


23 


18 


13 


12.47 


-17 54 28.26 


19.01 


15.37 


13.95 


10.65 


9.83 


9.44 


9.19 


9.21 


9.10 


9.14 


ukirt 


09-B3 


24 


18 


13 


23.87 


-17 53 10.10 








11.68 


10.75 


10.31 


9.96 


9.89 


9.80 


9.91 


ukirt 


09-B3 


25 


18 


13 


25.90 


-17 53 57.45 






20.28 


12.33 


11.34 


10.82 


10.48 


10.45 


10.25 


10.46 


ukirt 


09-B3 


26 


18 


13 


25.23 


-17 51 35.64 


15.37 


14.90 




13.64 


11.74 


10.84 


10.16 


10.12 


10.00 


10.07 


ukirt 


Q9-B3 



Note. — For each star, number designations and coordinates (J2000) are followed by magnitudes measured in different bands. J,H, and K B 
measurements are from 2MASS, while the magnitudes at 3.6 £tm, 4.5 /im, 5.8 ^tm, and 8.0 fim are from GLIMPSE. B, V, and _R associations 
arc taken from the astrometric catalog NOMAD. 

a 2MASS magnitudes arc above saturation limits, and have uncertainties ~ 0.3 mag 
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Fig. 2. — K-band spectra of early-type stars obtained with NIRSPEC on the Keck telescope. The 
full line shows the spectrum after correction for atmospheric absorption and instrumental response, 
and interstellar extinction. Extinction correction wa s performed with the values listed in Table EJ 
and by assuming a power law with an index of —1.9 (IMessineo et al.ll2005l ). The dashed line shows 
a black body with a temperature equal to the stellar effective temperature. Each spectrum was 
normalized at 2.12 /im, and arbitrarily shifted for clarity. The locations of detected spectral lines 
are indicated with marks and labels. 
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Fig. 3. — K-band spectra of early-type stars obtained with UIST on the UKIRT telescope. Curves 
and labels are as in Fig. [2] 
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Fig. 4. — Spectra of star $4, which is an already known WR star (jHadfield et al.ll2007l ). and star 
#7, which is a newly dete cted WR star. Fo r comparison, a spectrum of a WN7 star, WR120, is 
presented from the atlas of lFiger et al.l (|1997l ). The rest wavelengths of Hell, Hel, HI, and CIII/NIII 
lines are also marked. 
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and are likely to be members. The high degree of completeness of the spectroscopic observations 
of the bright sample allows for a precise numbering of stars in various post-main sequence phases 
(WR, cLBV, RSG, and OB stars). 



3.2.2. Color- Magnitude Diagram and cluster reddening 

In the top panels of Fig. [5j we display a (J—K s ) vs. K s diagram of 2MASS data points within 
the cluster radius of 3.5', as well as a diagram of a comparison field. The same CMD with deeper 
SofI data, covering the central 4.9' x 4.9', is shown in the bottom panel of Fig. [5J 

Several distinct populations of stars are seen in the CMD. There is a sequence of stars at 
J— K s =0.5 mag, which we attribute to field stars in the closer Sagittarius-Carina spiral arm. There 
is a broad sequence with J—K s from 3 to 6 mag, which is populated mainly by field giant stars in 
the background of the cluster. Three candidate RSGs were detected from this sequence. Finally, 
there is a sequence of bright stars with J— iT s =~ 1.5 mag and with K s from 3 to 14 mag (bottom 
panel of Fig. [5]). This sequence is seen only in the cluster field, and this is the brightest extension 
of the cluster sequence. 

To confirm and isolate the cluster sequence, we performed a statistical decontamination. 
2MASS data were used in order to have a comparison field. We counted the number of stars 
in a grid of 0.5 mag in (J—K s ) and of 1.0 mag in K s . From each bin of the cluster CMD, we 
subtracted a number of stars equal to that of the corresponding field bin. First, we subtracted all 
contaminating giant stars (spectroscopically detected), then we continued with a random subtrac- 
tion. The decontaminated diagram is displayed in the right panel of Fig. [U and different symbols 
indicate known spectral types. The cluster sequence appears populated by massive stars. An av- 
erage interstellar extinction of ^4k s = 0.83 ± 0.2 mag (Av= 9.1 mag) was estimated by matching 
the colors of the o bserved cluster sequence w ith a theoretical isochrone of solar composition from 



oc A 19 (IMessineo et al J 120051 ). This measurement is independent of age since the isochrones are 



the Gene va group (ILeieune fc Schaererll200lh . and by assuming a power-law extinction curve Ay 



almost vertical sequences in this plane. An isochrone of 4.5 Myr and solar me tallicity, taken from 



the n on-rotating models of the Geneva group, is over-plotted on the diagram (jLejeune &: Schaerer 



200l|). 



4. Luminosities of the massive stars 



In the following, we analyze the properties (colors, magnitudes, and luminosities) of massive 
stars in the CI 1813-178 cluster, which we list in Table [2j A kinematic distance of 4.8 kpc is used 
(see below the RSG section). For each star, we assumed an intrinsic c olor consistent with its stellar 
sp ectral type (IKoornneeJ Il983l ; iMartins h, Plea 120061 ; IWegnerl 1 19941 ) ; we used the extinction law 
by IMessineo et all (|2005l ) . and the kinematic distance, and we estimated a value of ^4k s and Mk 
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Fig. 5. — Left upper panel: A 2MASS J—K s versus K s diagram of point sources detected within 
3.5' from the cluster center is displayed. Middle upper panel: For comparison, A CMD for stars 
taken from a surrounding annulus (with a minor radius of 7' and area equal to that of the cluster) is 
show. Right upper panel: A CMD of the cluster field after statistical decontamination is shown. 
The vertical lin e indicates an isochrone o f 4.5 Myr, with solar metallicity, and with a distance of 
4.8 kpc (solid) (jLeieune fc Schaererll200ll ); the same isochrone shifted to a reddening of ^4k s =0.83 
mag is indicated with a dashed line. The locations of OB stars on the ZAMS with ^4k s =0 mag 
are marked with labels; dashed lines show their paths to ^4x8=0.83 mag. The arrow indic ates the 
reddening vector for ^4^=1 mag, based on the extinction law by iMessineo et al.1 (120051 ). Stars 
spectroscopically observed are marked with special symbols, as listed in the legend. Bottom: 
J—K s versus K s diagram of point sources detected in the SofI images. Symbols are the same as in 
the upper panels. 
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Fig. 6. — Upper panel: K s — [8.0] versus J—K s diagram of point sources detected within 3.5' 
from the cluster center. The arrow indicates the reddening vector for ^4k = 1-0 mag, which is 
calculated by following a near-infrared power law with a =1.9 (Mess ineo et al.l 120051). an d by using 
the extinction ratios A^.q/ Ak s , ^U.s/^Ksj A^^/Ak s , Ah nM^of llndebetouw et al.l (|2005l ). Symbols 
are as seen in Fig. Bottom panel: GLIMPSE [3.6]-[8.0] versus [3.6]-[5.8] diagram of point 
sources detected within 3.5' from the cluster center. Symbols and arrow are as shown in the top 
panel. 
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(see Table [2]). The Mk values were transformed into bolometric magnitudes, Mbol, by adding the 
bolometric corrections. The adopted effective temperatures and bolometric corrections per spectral 
type are given in the Appendix. 



4.1. Red supergiants 



Star #1 is the brightest star within the cluster area, K s =3.8 mag, and is consistent with being 
a K2-K5 I cluster member (Paper !) . From the radia l velocity (Vlsr = 62 ±4) of star $4 (Paperl), 
and the Galactic rotation curve bv lReid et afllboOflh. we obtained a revised near kinematic distance 
of 4. 8ln'7g kpc- With intr insic colors from Koornneei (1983), bolometric correction, B Ck, taken 
from iLevesque et al.l (|2005l ). we estimated a luminosity L^= 9 ± 7 x 10 L@. 



Stars #32, #38, and #39 have EW(CO)s larger than the typical values for red giants. Stars 
#32 and #38 were observed with NIRSPEC, and their i^-band spectral coverage allows us to study 
both the CO band heads, and the shape of the continuum. The large values of EW(CO) s and the 
absen ce of water absorption suggest a supergiant luminosity class for stars #32 and #38 (jComeron 
20041 ) . Star #39 was observed with UIST on UKIRT, with the K— long filter only, and we do not 
have information on the shape of its continuum. We derived spectral types M2.5I, M3.5I, and Mil, 
respectively, by assuming a supergiant luminosity class for all three stars. 

Stars #32, #38, and #39 are among the brightest stars with red J—K s (from 3 to 9 mag), 
much redder than the cluster sequence as shown in Figs. and [H By assuming intrinsic H—K s 
and J—K s colors for a given spectral type (see the Appendix), we obtained Ak s = 2.3 mag, 3.6 mag, 
and 2.3 mag, respectively. These values of interst ellar extinction are th e combined contribution of 
interstellar and possible circumstellar extinction (jMessineo et al.ll2005l ). By analyzing the distri- 
bution of detected stars in the Ks — [8.0] vs J—K s plane ©, star #38 appears to show an excess 
in Ks — [8.0], which is suggestive of a circumstellar envelope. Stars #32, and #39 have Ks — [8.01 



and J—Ks colors in agreement within errors with those of naked late-type stars (jMessineo et al 



2004 



20051 ). We corrected for t he total extinct i on by using the intrinsic colors for given spectral 



types and the extinction law by lMessineo et al.l (|2005l ). We then artificially reddened the stars to 
emulate the average cluster condition. By assuming these stars were moved along the reddening 
vector to the average cluster extinction (^4x5=0.8 mag), we obtained K s =6.7 mag, 5.8 mag, and 
7.2 mag, respectively. 

If these RSGs were at the same distance as the stellar cluster they should be the brightest 
stars at near-infrared wavelengths, and their K s magnitudes would increase with later spectral 
types. Their spectral types are later than the bright (K s =3.8 mag) K2I star, and their de-reddened 
magnitudes are similar to those of the brightest OB stars (K s =~ 7 mag). The spectral types and 
magnitudes of these three RSGs are not consistent with cluster membership; they are too faint to 
be RSGs and part of the cluster. This reasoning is valid also for stars with moderate mass-loss, 
i.e. #38, because bolometric corrections BCk are almost constant for colors -Kg— [8.0] < 3 mag 
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(chapter5, Figs. Al and A2. iMessineo 112004 ). 

We estimated a density of non-member RSGs equal to 0.16 in the central 3.5'. These RSGs 
could represent an outburst farther away, perhaps related to the Galactic bar. 



4.2. Wolf Rayet stars 



Star #4, a WN7b, was spectroscopically detected by lHadfield et al.l (|2007l ). Hadfield et al. 
performed a successful color-based selection of candidate WR stars with 2MASS and GLIMPSE 
data. WRs have an infrared excess (see Figs. [6] and [5]), which is caused by free- free emission, and 
do not follow the reddening vector. Star #7 is a newly detected WN7o star, its colors fall within 
the Hadfield et al. criteria. 

For stars #4 and #7, we used int r insic c olors for WN stars with broad and narrow lines, 
respectively, derived by ICrowther et al.l (|2006al ) (see the Appendix) . We obtained an interstellar 
exti nction Ak„=0.7 mag, an d ^4k s =0.8 mag, respectively. By adopting the kinematic distance, BCk 
from lCrowther et al.l (j2006ah . we obtained L*= 5.8±2.5xl0 5 L Q for star #4, and U= 4.8±2.1xl0 5 
L & for star #7. 



4.3. Luminous blue variables 



Among the observed spectra, we detected one candidate cLBV s tar (#15). LBV stars are rare 



mass i ve supergiants in transition t owards the Wolf-Rayet phase (e.g. IConti et al. 



1995 



Nota et al 



1995 



Figer et al.l[l997l : IContilll984h . Since their evolutionary paths on the CMD are quite uncertain, 



their detections in stellar clusters are of primary importance. The K-band spectrum of star $45 is 
similar to that of a P Cygni-type B supergiant, cLBV, e.g. the [OMN2000]LS1 star presented by 



Clark et al 



(2009j). We u sed CMFGEN, the iterative non-LTE line blanketing method presented 



bv lHillier Milled (119981). to estimate the physica l pro perties of this super giant. For details on the 
modeling see lNaiarro et al.l (|1999l ). iNaiarrol (|200ll ) and Iciark et al.l (|2009l ). 



The K-band spectra provided the primary diagnostics (the Hel lines at 2.058 /mi and 2.112 
/xm, and the Mgll lines at 2.138 /iin and 2.144 fim). We obtained an effective temperature T e jj =~ 
16, 000 ± 300 k, and a value of L*= 2.1 ± 1.5 x 10 5 Lq. This value matches the luminosities of the 



faintest P-Cygni supergiants known in the Milky Way (jClark et al.ll2005l . |2009| ). 



4.4. OB stars 



We detected 21 OB stars in the CI 1813-178 cluster. The current dataset of medium-resolution 
K-band spectra allow for a spectral classification of OB stars typically within ±2 spectral types 
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(jHanson et al.lll996l ). The combination of spectral and photometric properties suggests that we 
have detected 14 OB supergiants, 4 OB giants, and 3 OB main sequence stars. 

For star #6 (06-0 7If+), when using an intrinsic (H—K s )q = —0.1 mag, and BC^=— 4.01 mag 
(|Martins fe Plea 1 hood ) . we obtained M K =-5.87 mag and L*= 8.2 ± 2 .3 x 10 5 L r .,. The estimated 



value of Mu is co nsistent with the value given for an 06-07 star by iMartins &: Plea ([2006) and 
Clark et al.l l|2005l ). 



Star #16 is anoth er rare O supergiant, a n 08-09K type. By assuming (H—K s )q = —0.1 mag, 
and 5C K =-3.84 mag (IMartins fc PlezlboOfih . we obtained M K =-7.32 mag and L*= 23 ± 6 x 10 5 
L & . Star #16 is the brightest early-type cluster member, and it is loca ted close to the o bserved 
Humphreys-Davidson limit for stars with similar effective temperatures (IClark et al.1 12005). 



From a comparison with the Mk values by IMartins &: Plea (|2006l ) and iPanagial (j 19731 ) , the OB 
stars #25, #24, and #26 are likely dwarfs, since they have Mk< —4 mag (£*< 1.4 x 10 5 L Q ). The 
OB stars #9, #21, #22, and #23 are probably of luminosity classes III or II, since they are fainter 



than a typical BSG. Their Mk values range from —4.7 to 
1.5 x 10 5 Lq). All remaining OB stars are supergiants. 



-5.36 mag (L* varies from 1.1 x 10 5 to 



The large spread in magnitudes of OB supergiants (~ 3 mag in if s -band) is not surprising. A 
similar observatio nal spread is observed in the 2MASS K s magnitudes of BSGs in Westerlund 1 
ddark et alJbood ). 



4.5. X-ray emitters 



X-ray emission m ay be generated in th e circumstellar envelopes of massive stars due to their 



strong shocked winds (iLucv & : White Il980l ). Single OB stars emit with a typical X-ray luminosity 
of Lx= 10 31-33 erg s _1 ( Pollocklll987l ). and have a typical ratio between the X-ray and bolometric 
luminosities of about 10~ 7 . 

1Q 32- 



OB+OB a nd OB+WR binaries generally have higher luminosities 



(Lx= 10 erg s . IClark et al.l 12008) . X-ray emission enables us to characterize the physical 
condition of stellar atmospheres, and to identify binary systems. 

X-ray observations of the CI 1813-178 cluster region were performed by Funk et al. (j2007i ) 
and lHelfand et al.l (120071 ). successfully detecting a large number (75) of X-ray emitters. We looked 
for possible associations between X-ray emitters and massive members of the CI 1813-178 cluster 
(Paperl). All but one X-ray sources with a bright 2MASS counterpart (K s < 9.6 mag), were found 
associated with early- type cluster members (see Table [3j). Two X-ray emitters coincide with WR 
stars (#4 and #7); six others are associated with OB stars (#2, #3, #5, #6, #8, and #9). The 
remaining X-ray emitter (#71. lHelfand et al.ll2007l ) coincides with star #10, which is likely a cluster 
non-member. This star is located 6.6' from the cluster center, outside of the cluster radius (3.5'). 
Its spectrum has CO bands at 2.29 /im, which indicate a late-type star. 



We estimated the ratios between the X-ray and bolometric luminosities, and we compared 



Table 2. Derived properties of stellar members. 



ID 


Ks 


J-Ks 


H-Ks 


Ej- 


Ks 


Eh-Ks 


A K 




Lum 


log(Teff[K]) 


Sp. type 


Class 


1 


3.79 


1.67 


0.44 


1 


02 


0.31 


0.54± 


0.13 


-10.16 


4.95 


24 
-0.25 


3.60 


01 
-0.02 


RSG 


1 


2 
3 


7.22 
7.79 


1.42 
1.41 


0.50 
0.46 


1 
1 


50 
49 


0.54 
0.50 


0.80± 
0.79± 


0.02 
0.02 


-6.99 
-6.41 


5.70 
5.46 


34 
-0.34 

-0.34 


4.31 
4.31 


0.12 
-0.12 
12 
-0.12 


B0B3 
B0B3 


1 
1 


4 


7.94 


1.68 


0.66 


1 


31 


0.39 


0.69± 


0.05 


-6.16 


5.76 


n ifi 
-0.17 


4.70 


-0.05 


WN7 


1 


5 


8.56 


1.40 


0.50 


1 


61 


0.60 


0.87± 


0.02 


-5.72 


5.75 


12 
-0.13 


4.50 


03 
-0.03 


0709 


1 


6 


8.57 


1.72 


0.58 


1 


93 


0.68 


1.03± 


0.02 


-5.87 


5.92 


11 
-0.12 


4.54 


02 
-0.02 


0607If 


1 


7 


8.66 


1.64 


0.61 


1 


51 


0.50 


0.80± 


0.05 


-5.55 


5.68 


16 
-0.17 


4.70 


04 
-0.05 


WN7 


1 


8 


8.75 


1.11 


0.37 


1 


32 


0.47 


0.71± 


0.01 


-5.36 


5.61 


0.12 
-0.13 


4.50 


0.03 
-0.03 


0709 


1 


9 


9.34 


1.40 


0.46 


1 


61 


0.56 


0.86± 


0.02 


-4.93 


5.47 


13 
-0.14 


4.51 


03 
-0.03 


0709 


III 


11 


9.59 


2.15 


1.26 


2 


23 


1.30 


1.37± 


0.02 


-5.19 


4.98 


-0.34 


4.31 


12 
-0.12 


B0B3 


III 


12 


6.75 


2.32 


0.75 


2 


23 


0.74 


1.19± 


0.02 


-7.85 


5.07 


n 

u. io 
-0.14 


3.98 


u.uo 
-0.03 


B9A2 


1 


13 


6.85 


1.94 


0.71 


2 


02 


0.75 


1.08± 


0.03 


-7.64 


5.96 


U.04: 

-0.34 


4.31 


12 
-0.12 


B0B3 


1 


14 


6.91 


1.30 


0.47 


1 


38 


0.51 


0.74± 


0.02 


-7.24 


5.79 


0.34 
-0.34 


4.31 


0.12 
-0.12 


B0B3 


1 


15 


6.96 


1.97 


0.82 


1 


99 


0.84 


1.09± 


0.02 


-7.53 


5.38 


22 
-0.23 


4.20 


01 
-0.01 


LBV 


\ 


16 


7.33 


2.10 


0.77 


2 


31 


0.87 


1.24± 


0.01 


-7.32 


6.36 


nil 

-0.12 


4.48 


n ni 

U.U-L 

-0.01 


0809If 




17 


8.52 


1.18 


0.45 


1 


26 


0.49 


0.68± 


0.03 


-5.56 


5.13 


34 
-0.34 


4.31 


12 
-0.12 


B0B3 


1 


18 
19 


8.61 
8.72 


1.44 
1.53 


0.44 
0.53 


1 
1 


65 
74 


0.54 
0.63 


0.88± 
0.94± 


0.05 
0.02 


-5.67 
-5.62 


5.74 
5.72 


0.13 
-0.14 
0.12 
-0.13 


4.50 
4.50 


0.03 
-0.03 
0.03 
-0.03 


0709 
0709 


1 


20 


9.14 


1.46 


0.47 


1 


62 


0.51 


0.86± 


0.03 


-5.13 


5.58 


0.11 
-0.12 


4.54 


0.01 
-0.01 


0708 




21 


9.29 


1.44 


0.49 


1 


60 


0.57 


0.86± 


0.03 


-4.98 


5.04 


0.38 
-0.38 


4.36 


0.13 
-0.13 


09B3 


III 


22 


9.34 


2.08 


0.76 


2 


24 


0.84 


1.22± 


0.03 


-5.28 


5.16 


0.38 
-0.39 


4.36 


0.13 
-0.13 


09B3 


III 


23 
24 


9.44 
10.31 


1.21 
1.37 


0.39 
0.44 


1 
1 


34 
50 


0.42 
0.47 


0.71± 
0.79± 


0.02 
0.02 


-4.67 
-3.89 


5.03 
4.72 


0.30 
-0.30 
0.30 
-0.30 


4.39 
4.39 


0.11 
-0.11 
0.11 
-0.11 


09B3 
09B3 


III 

V 


25 


10.82 


1.51 


0.52 


1 


64 


0.55 


0.88± 


0.04 


-3.46 


4.55 


0.30 
-0.30 


4.39 


0.11 
-0.11 


09B3 


V 


26 


10.84 


2.80 


0.90 


2 


93 


0.93 


1.55± 


0.02 


-4.12 


4.81 


0.30 
-0.30 


4.39 


0.11 
-0.11 


09B3 


V 
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Table 3. Massive stars with X-ray emission. 



ID 


K s 


Sp. Type 




Countsc/i 


HR 


UxMM 


CountsxMM 


Lx 




[mag] 














[10 31 erg s" 1 ] 


2 


7.22 


OB 


39 


14.10 


-0.33 






4.1 


3 


7.79 


OB 


43 


16.40 


-0.53 






4.8 


4 


7.94 


WR 


24 


272.40 


0.78 


2 


238 


80.5 


5 


8.56 


OB 


41 


214.50 


0.42 


4 


138 


63.4 


6 


8.57 


OB 


58 


13.60 


-0.57 






4.0 


7 


8.66 


WR 


37 


34.90 


0.89 






10.3 


8 


8.75 


OB 


27 


13.00 


-0.82 






3.8 


9 


9.34 


OB 


36 


9.50 


-0.20 






2.8 


10 


9.60 


K 


71 


71.10 


-0.31 






21.0 



Note. — For each star, number designations, K s magnitudes, and spectral types from Tables Q] 
and [7] are followed by the nu mber designat i ons (I Dch), counts (counts), and hardness (HR) from 
the Chandra observations by Helfan d et al.l (|2007l ) . and by the number designations (Hxmm) arid 
XMM counts reported by Funk et al.1 ( 2007 ). Estimates of Lx are taken from Paperl, and assume 
a distance of 4.7 kpc, N(H)=1.6 x 10 22 cm -2 , a power-law model, and a photon index of 1.5. X-ray 
luminosities are given in units of 10 31 erg s _1 . 
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the ratio, hardness and X-ray luminosities of our nine X-ray emitters with those of Chandra point 
sources detected in Westerlund 1 (Fig. 5 IClark et al.1 120051 ). In Westerlund 1, WR stars have 
luminosities Lx larger than ]_o 32-34 erg s _1 and hardness from —0.1 to 1, while most of the detected 
OB stars have Lx=10 32 erg s _1 (which is consistent with a ratio of ~ 10~ 7 ) and hardness between 
—0.8 and —0.1, as expected for single stars with shocked winds. Clark et al. suggest that all 
OB stars with X-ray emission significantly harder than —0.5 in Westerlund 1 are binary systems. 
This conclusio n is supported b y emerging ey i dence of a high-binary fraction of massive stars in 
Westerlund 1 (jBonanos 1 120071 ). ISkinner et al.1 (|2010l ) propose a number of other possible scenarios 
to explain the X-ray emission of single WN stars. Magnetic wind confinement could also explain the 
presence of a hot plasma component without invoking the presence of a close companion. However, 
current detections of magnetic fields in WN stars are absent. 

Star #4, a WN7 star , is the brightest X-ray emitter. It is associated with the Ch a ndra source 
#24 (jHelfand et al.ll2007h . and coincides also with the XMM source #2 of iFunk et all (j2007h . The 
ratio between the X-ray and bolometric luminosities of star #4 is 3.7 x 10~ 7 . For star #7, we 
obtained a ratio 7 times fainter. Both ratios are consistent with those measured in other WN stars 



by ISkinner et al.l (|2010l ). Their high values of hardne ss (0.78 and 0.89 in Table [3|) are consistent 
with those of colliding w ind binaries (IClark et alj|2008r ). Since only a few other late WN stars have 
been detected in X-ray (jSkinner et al.ll2010i ). our new detections are a significant addition. 



Star #5 (07-09) was detected by both the XMM and Chandra satellites (see Table EJ). It 
is a strong X-ray emitter (Lx= 1 x 10 32 erg s _1 ); the ratio between the X-ray and bolometric 
luminosities is 2.9 x 10~ 7 . Besides the two WR stars, star #5 is the only other source with a strong 
X-ray hardness (0.4). The high values of X-ray luminosity and hardness indicate that #5 is another 
binary. 

For the BSGs #2 and #3 (B0-B3), we measured ratios between the X-ray and bolometric 
luminosities of 2.1 — 4.3 x 10~ 8 , which are in agreement with the ratios measured for single stars 
later than Bl (ICohenl Il99fil : Iwaldron fc Cassinellil 120071 ). Stars #8 and #9 (07-09) have also a 
ratio of about 2.5 x 10~ 8 . 



For star #6, an 0607If, we estimated a ratio of 1.3 
(—0.57) suggest radiative shocks in stellar winds. 



10 8 . The low ratio and low hardness 



5. Spectro-photometric distances 



Massive evolved stars in a young stellar cluster span a range of masses (see Fig. [7|), therefore, 
of luminosities. Spectro-photometric distances from OB supergiants are less accurate than those 
from dwarfs an d giants. When considering the dwarfs and the My; magnitudes fo r 09, BO, Bl, 



and B2 dwarfs (Martins fc Ple2 



2006; 



Humphreys McElrov 



1984: 



Koornneeall983l ). we estimated 



distances of 3.8 ± 0.6 kpc, 3.2 ± 0.5 kpc, 2.6 ± 0.4 kpc, and 1.9 ± 0.3, respectively. 
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From the photometric prop erties of the candidate g i ants (#9, #21, #22, and #23), and Mk 
values for classes II and III from Humphreys fc McElrovl ()1984l ) and IWegnerl (j 19941 ). we derived an 
average spectro-photometric distance of 3.0 ± 0.7 kpc for class III, or 4.5 ± 0.6 for class II. 



For the two WN 7 stars, we adopted in trinsic magnitudes from ICrowther et al.l (|2006al ) , and 
the extinction law by lMessineo et al.l (120051 ) . and obtained a distance of 2.6 kpc for star #4, and 
5.8 kpc for star #7. The average distance for the two WN7 stars is 4.2 ± 1.6 kpc. 

The derived spectrophotometric distances are listed in Tabled While the distance estimates 
from giants and WRs within errors are consistent with the kinematic distance, the distance esti- 
mates from dwarfs are only consistent with the kinematic distance if the dwarfs are all late O stars. 
Higher-resolution spectra are needed to refine the spectral types. 

Observa tions of radio hydrogen recombinatio n lines of the W33 complex reveal two velocity 



1978 



components (|Bieging et al 
kms -1 , and the rotation curve of 
3.53 



Goss et al 



Reid et al 



-O.40 

62 kms -1 



-o 46 kp c > while we obtained a distance of 4.82_ 2g 



1978|) . By using the radial velocity component at 35 
(2009), we calculated a near kinematic distance of 
~ rU " 1 kpc with the radial velocity component at 



A radio monitoring program of methanol masers in the direction of the W33 complex is 
currently ongoing. It will yield parallactic distances of the masers (Brunthaler et al. in preparation). 



6. Progenitor masses and the cluster age 

The CI 1813-178 cluster contains a large number of evolved massive stars (BSGs, a cLBV, 
two WN7 stars, and one RSG star), which are listed in Table [TJ We plot the inferred stellar 
luminosities versus stella r effective temperatures in Fig. [HJ together with theoretical stellar models 



from the Geneva group (jMevnet &: Maederll2000l ). By comparison with the models, we estimated 



initial masses from 20 to 100 M@. 

Models predict a large span of stellar masses in post main sequence phase, e.g. a population 
of 4.0 Myr would have a mass of 25-35 Mq at a TO, but it would still contain stars of about 50-M©, 
or even 100 if rotating models are considered. This is well illustrated in Figure where we show 

Table 4. Average spectrophotometric distances 



Sp. type distance 



OBV 2.9 ±0.8 
OBIII/II 3.8 ±1.0 
WN7 4.2 ± 1.6 



average 



3.6 ±0.7 
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Fig. 7. — For a given age, the mass at the TO and the maximum m ass predicted by models 
are shown. Several models are used: a model b y ISchaller et al.l (119921 ) with a solar metallicity 
(asterisks), model s with metallic i ty tw ice solar by ISchaerer et al.l (j 19931 ) (diam onds), models with 
high m ass-loss by lMeynet et al.l (119941 ) (triangles), and the rotating models by lMeynet &: Maeder 
(J2000) (crosses). For each model, the maximum masses are connected with a dotted line, while the 
masses at the TO are connected with a dashed line. 
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the predicted maxi mum initial ma s ses an d T O masses as a function of cluster age. We used non- 
rot a^iiig_jmojiels by]Sd (|1992l ) and ISchaerer et al.l (11993T). models with high-mass loss 
by iMevnet et al.l (|1994l ) , and the rotating models by iMevnet Maederl (|2000l ) . Metallicity ranges 
from solar to twice solar. For the same age, different models predict variations in the initial masses 
of up to 50%. 

The mixture of evolved stars indicates progenitor masses larger than 20 M@, and more likely 
of 25-60 Mq, with a few exceptions. 

The Of stars are among the mo st luminous stars. The lu minosity of the 08-09If star suggests 
an extremely massive star (100 My. lMeynet fc Maederl l2003). The 06-07If star has an estimated 
U= 8.1 ± 2.2 x 10 5 Lq, which is predicted for a 30-70 Mq. 

WR stars of WN7 type have been found only in stellar clusters with masses at the TO larger 



than 35-40 M Q ( Magsey. et al. 



2001 



Clark et al.ll2005l ). The lumimosities of the WR stars in the 
CI 1813-178 cluster (L*= 5.2 ± 2.3 x 10 5 and 4.5 ± 1.9 x 10 5 Lq) are similar to those inferred for 
WN7 stars in the Westerlund 1 cluster (L* from 3.1 x 10 5 to 5.0 x 10 5 Lq). The stellar luminosities 
suggest initial masses from 40 to 70 Mq. 

The cLBV and RSG have luminosities expected for less massive stars (~ 20 Mq). However, 
the 2MASS magnitudes of the RSG star have errors of about 0.3 mag, due to saturation, and the 
resulting luminosity could be underestimated. 



Models of young simple stellar populations by lMeynet &; Maederl (|2003l ) predict that stars with 
masses between 9 to 25-35M have a RSG phase. Single WR stars have initial mass es greater than 
26 — 30 Mq, while binary WR stars have initial masses greater than 20 — 25 Mq (jEldridge et al 



20081 ). By assuming coevality between the RSG and the WR stars, the CI 1813-178 cluster would 
be between 4 and 6 Myr years old. A cluster with stars exceeding 100 Mq (like the 08-09If) 
would require an age of 3-4 Myr, and a turn-off (TO) at 35 Mq. Therefore, the simultaneous 
presence of RSGs, WRs, and Of stars, would further narrow the possible age range to 4-4.5 Myrs. 
However, either the luminosities of the cLBV, of the RSG star, and/or of the late B supergiant 
are underestimated by ~ 0.7 dex, or the 'cluster' is more of an 'association' with some degree of 
non-coevality. 

The cLBV star appears to be rather faint. The estimated lum inosity of 2.4 x 10 5 Lq is similar 
to that of the HD168607 and HD316285 LBVs (jciark et alJhood ). Other known LBVs in clusters 
are typically among the m ost luminous members, e.g. qF362 and the pistol star in Quintuplet 
(e.g. iMauerhan et al.ll2010l ). In Westerlund 1 the W243 LBV has an initial mass of about 40 Mq, 
which is consistently similar to the masses of the cluster WRs (40-50 Mq), and larger than the 
estimated mass at the TO (30 M Q ) (jRitchid \20(M ). The cLBV in C11813-178 has a mass smaller 
than that of the two detected WR stars. Its mass is consistent with the mass of the RSG, and 
would require a mass at the TO of about 20 Msun, and an age of 5-7 Myr. Further observations are 
recommended to explore the degree of coevality in C11813-178. Some non-coevality would easily 
explain the discrepant masses. A population with an age of 4.0 — 4.5 Myr, and a spread in age of 
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1 Myr, could explain the observed range of stellar masses. 

The stellar mass at the TO for a population with an age of 4.0 — 4.5 Myr is likely between 
25-35 Mq. This is consistent with the observations of dwarfs and giants. Three OB dwarfs with 
K s < 10.0 mag were detected with 09-B3 spectral types, while the giants with 07-B3 types have 
if s =9.3-9.5 mag. An 07V star (~ 32 Mq) at an a verage extinction of ^4 k„=0-8 mag, and a distance 
of 4.8 kpc is expected to have a K s =~ 10.4 mag ([Martins fc PlezJhood ). while an 09V (~ 22 Mq) 



star has a if s =10.9 mag. 

A further spectroscopic survey of fainter stars is needed to sample the TO region. Moreover, 
high-resolution spectra would allow for a more precise spectral classification. By narrowing the 
errors shown in Fig. El and increasing the sample, we will be able to better constrain the age, and 
verify the degree of coevality. 



7. The cluster mass 



By adding the masses of the spectroscopically identified massive stars, we estimated a minimum 
cluster mass of 990 Mq. We considered all 34 stars with J—K s between 1 and 3 mag, and K s < 10 
mag. By assuming that they all have masses greater than 35 Mq, and using a Salpeter mass function 
integrated from 0.8 to 120 Mq, we estimated a cluster mass of 13000 ± 3000 Mq, where the error 
is the Poisson error of the number of massive stars. By assuming that these 34 stars have masses 
greater than 25 Mq, we obtained a cluster mass of 8700 ± 2000 Mq. The two calculations take 
into account the uncertainties of the mass at the TO (25-35 Mq), predicted for a coeval population 
with an age of 4-4.5 Myr. 

The CI 1813-178 cluster is, therefo re, a new addition to the list of 13 known y oung massive 



clusters (> 10 4 Mq) in the Milky Way (jMessineo et al 



2009 



Negueruela et al.ll2010l ). 



8. Cluster surroundings 



A 24 /im image from MIPSGAL, the survey of the inner Galactic plane using the M ultiband 
Infrared Photometer for Spitzer aboard the Spitzer Space Telescope (jCarey et al.ll2009l ). of the 
whole W33 region is shown in Fig. [9l together with contours of radio continu um emission at 20 cm 
(jWhite et al.ll2005l ). The complex extends over an area of roughly 25'x20' (jBieging et al.lll978l ). 
Radio observations show that W33 is made up of a number of d i screte sour ces. Some of these ra dio 
sources have been classified as candidate SNRs bv lBrogan et al.l (12006T) and Helfan d et al.l (|2006l ) on 
the basis of morphology and spectral indexes, with MAGPIS data ([White et al.ll2005l ). See Table 

El 

The CI 1813-178 cluster appears located on the edge of the W33 complex. The radial velocity 
of the K2I star in the CI 1813-178 cluster is of 62±4 kms" 1 (Paperl), and well agrees with the 
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Fig. 8. — Luminosities of massive stars in the CI 1813-178 cluster are plotted versus effective 
temperatures. Values are taken from Table [2j Stellar tracks for stars of 20, 25, 40, 60, and 85 Mq, 
based on rotating Geneva models with a solar metallicity, canonical mass-loss rates, and initial 
rotational velocity of 300 kms -1 are also shown with dotted, dashed and continuous lines. For 
clarity, we use different line styles together with labels to indicate each model. 
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Fig. 9. — 24 yum emission from MIPSGAL (jCarey et al.ll2009i ) of th e whole W33 com plex. The 
white contours show 20 cm radio continuum emission from MAGPIS (j White et al.ll2005l ). Contour 
levels are 0.003, 0.009 and 0.020 mJy beam" 1 ; the beam size is 6.2" x 5.4", FWHM, and the position 
angle of the major axis is along the north-south direction. Circles and labels indicate the location 
of the stellar clusters and candidate SNRs. 
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high-velocity gas component of W33. The spatial coincidence of the CI 1813-178 cluster with the 
SNR G12.82-0.02 and G12.72-0.00 has already been reported in Paperl, and support the association 
of the cluster with the complex. The filamentary shape of the CI 1813-178 stellar cluster and its 
location on the edge of W33 suggests a secondary episode of star formation, perhaps triggered by 
an expanding shell. 

We searched for other possible candidate clusters and associations with SNRs in the direction 
of the W33 complex. A number of candidate stellar clusters have been identified in GLIMPSE 



and 2 MASS images in the direction of the W33 reg ion by IMercer et al.l (|2005l ) and iBica et al . 
J200I), which we lis t in Table El The W 33 MYSO (jDavies et al.l bold ) is projected into SNR6 
(G13. 1875+0.0389) faelfand et al.l bood ). The spatial coincidence of the SNR and the MYSO 



suggests a physical association. The W33A MYSO could be an episode of triggered star formation 
in duced by a supernov a explosion. Mercerl candidate cluster is the object number $d in the list 
of IMercer et al.l (|2005l ). It was identified as a stellar overdensity in the GLIMPSE catalog with 



an automatic algorithm. It is located at the center of the molecular complex, and it appears as 
a spread overdensity. Two other candidates are reported in literature in the surrounding of the 
W33 complex, but without a clear connection with the W33 complex. The BDS2003-115 candidate 
cluster is about 20' North of the main W33 complex and is associated with SNR5 (G12.83-0.02) 
(jHelfand et al.ll2006l ). The BDS2003-7 cand idate cluster appears as a small group of stars without 
associated radio emission (IBica et al.l 120031 ) . 



We visually inspected the 2MASS images, and located two other clumps of stars (ell and cl2). 
The ell candidate appears as a group of point sources on bright nebular emission in the 2MASS 
K s image. Inspection of the GLIMPSE and MAGPIS images reveals the presence of an HII region, 
suggesting the presence of massive stars (see Tableland Figs. PTOl and [TTj) . The cl2 candidate 
is another small concentration of bright stars (K s =8-10 mag) in another HII region. Nothing is 
reported in previous literature about both, ell and cl2, clumps. In addition, we searched for stellar 
over-densities in the direction of the W33 complex using both 2MASS and GLIMPSE star-counts. 
Detections are hampered by strong variations of the background level, which are due to variations 
of interstellar extinction and nebular emission. A spectroscopic and photometric follow-up study 
of these regions with SINFONI and UKIDSS data is ongoing, in order to confirm the presence of 
massive stars. 

Near-infrared spectroscopic follow-up observations are needed to characterize these sources, 
and to confirm their association with the W33 complex. However, the associations of CI 1813-178, 
ell, and cl2 with HII regions and/or SNRs suggest that these are other condensation of massive 
stars in the W33 complex. 
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Table 5. List of candidate clusters in the direction of the W33 complex. 



ID 


RA 


DEC 


Radius(') 


References 


CI 1813-178 


18 13 24 


-17 53 31 


3.5 


Paper I 


Mercer 1 


18 13 57 


-17 56 46 


2.3 


(Mercer et al. 20051 


BDS2003-115 


18 14 05 


-17 28 29 


1.2 


CBica et al. 2003 


BDS2003-7 


18 10 55 


-17 41 25 


0.6 


CBica et al. 2003 


ell 


18 14 22 


-17 56 10 


1.8 


Present work 


cl2 


18 14 08 


-18 00 15 


1.0 


Present work 




Arc Minutes 

Center: R.A. 18 14 05.00 Dec -17 28 29.0 



»* ■•••*». i-» 

• •••• *.- •. . •■ 



Arc Minutes 

Center: R.A. 18 10 55.00 Dec -17 41 25.0 




Arc Minutes 

Center: R.A. 18 14 22.00 Dec -17 56 10.0 




2 1 0-1-2 

Arc Minutes 

Center: R.A. 18 14 08.00 Dec -18 00 15.0 




2 1 0-1-2 

Arc Minutes 

Center: R.A. 18 13 57.00 Dec -17 56 46.0 



Fig. 10.— 2MASS K s images of the BDS2003-115 cluster (top-left), the BDS2003-7 cluster (top- 
middle), the ell cluster (top-right), the cl2 cluster (bottom-left), and the Mercer 1 cluster (bottom- 
middle). 
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Fig. 11 — Composite images of SNR5, SNR6, ell, cl2, and CI 1813-178. In blue the GLIMPSE 3.6 
/xm, in green the 5.8/xm, and in red the 8.0nm. Ove r-plotted are contours of 90 cm radio continuum 
emission detected by MAGPIS (jwhite et alJbood ). 



Table 6. List of candidate SNRs 



ID 


Name 


RA 


DEC 


Radius (') 


References 


SNR1 


G12.26+0.30 


18:11:17 


-18:10:00 


4 


fBrogran et al. 2006: Helfand et al. 2006) 


SNR2 


G12. 58+0.22 ? 


18:12:14 


-17:55:00 


5 


fBrogran et al. 2006) 


SNR3 


G12. 72-0.00 


18:13:18 


-17:55:42 


4 


fBroean et al. 2006: Helfand et al. 2006) 


SNR4 


G12.83-0.02 


18:13:35 


-17:50:30 


2 


fBroean et al. 2006: Helfand et al. 2006) 


SNR5 


G13.1875+0.0389 


18:14:07 


-17:28:43 


2.5 


CHelfand et al. 2006) 


SNR6 


G12.9139-0.2806 


18:14:44 


-17:52:47 


1.5 


CHelfand et al. 2006) 
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9. Summary 

A near-infrared spectroscopic survey of the brightest stars in the direction of the CI 1813-178 
cluster is presented. Among the 61 observed stars, 25 massive stars were detected. Two WR stars 
of type WN7, a cLBV, and 21 OB stars were identified. Among the OB stars, a 08-09If star and 
a 06-07If star were discovered. Eight of these evolved stars also have X-ray emission, as detected 
by the Chandra and XMM satellites. The hardness of the X-ray emission from the two WN7 stars 
strongly suggests binary systems. 

A spectro-photometric analysis of the OB stars reveals 14 supergiants, 4 giants, and 3 dwarfs. 
From the giants, dwarfs and WRs, we derived average spectrophotometric distances of 3.8 ± 1.0 
kpc, and 2.9 ± 0.8 kpc, and 4.2 ± 1.6 kpc. The distances from giants and WRs is in agreement with 
the kinematic distance. The distance estimates from dwarfs are only consistent with the kinematic 
distance if the dwarfs are late O stars. 

The mixture of evolved massive stars is reminiscent of other Galactic young massive clusters, 
such as Westerlund 1, Quintuplet, Galactic center, and CI 1806-20. We estimated stellar luminosi- 
ties, therefore, masses by comparing the luminosities with evolutionary tracks from the Geneva 
group. By assuming a Salpeter mass function, we obtained a cluster mass of 1.0 ± 0.2 x 10 4 M Q . 
A likely cluster age of 4-4.5 Myr is derived, however, a spread in age of about 1 Myr cannot be 
excluded. In order to better constrain the degree of coevality, further spectroscopic observations 
are required. 

The CI 1813-178 cluster is located on the Western edge of the W33 complex. We have located 
several other candidate stellar clusters that could belong to the same complex. 
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A. Late-type stars 



The EW(CO)s can be used for classification in subclasses of late-type stars. There is a linear 
correlation between stellar effective temperatures and EW(CO)s of the CO bands. Furthermore, 



Fieer et al. 


2006: 


Davies et al. 


2007) 



Among the observed 60 stars, 37 are found to be late-type stars. For 36 of them, spectra were 
taken with the K-long grism, covering the region of the CO band-head at 2.29 /im. CO bands at 
2.29 //m in absorption are detected in all spectra. Star #55 was observed only with the K-short 
grism; however, the presence of Na lines (at 2.2075 and 2.2077 /tm) and Mg I at 2.11 /im suggests 
a spectral type later than GO. 

The equivalent widths of 32 of these late-type stars are compatible with that of giant stars 
with spectral type from K0III to M7III (Table [7]). For star #58 we do not report any spectral type, 
since its spectrum has a poor signal to noise. 

Star #32, #38, and #39 have EW(CO)s typical of RSG stars, with spectral types M2.5, M3.5 
and Ml, respectively. However, their photometric properties indicate that they are unrelated to 
the stellar cluster. Star #1 is consistent with being a K2-K5I cluster member (see Paperl). 

Young massive clusters with ages from 4 to 30 Myrs may contai n yellow supergia nts (YSG) 
stars, e.g. the Westerlund 1 (jciark et alJboOfJ ) and RSGC1 clusters (jFiger et alihood ). YSG are 
rare F or G-type supergiants in transition towards the RSG phase or, back from the RSG locus, 
evolving blue-ward. In a coeval population, YSGs or RSGs are expected to be brighter in K s than 
early-type massive stars. In the CI 1813-178 cluster, all detected stars with CO band-head at 2.29 
//m appear fainter than the brightest OB stars, and the gap between the brightest early-type star 
(i^ s =6.75 mag) and the K2I star (if s =3.79 mag) is devoid of stars. Furthermore, the CO band- 
heads of the five late- type stars with K s between 7 and 8 mag indicate late M giants (M2-M5), or 
late K supergiants (K2-K4). Their EW(CO)s are not consistent with being F or G supergiants. 



B. Effective temperatures and bolometric corrections 

In order to estimate stellar luminosities, estimates of effective temperatures and bolometric 
corrections as a function of spectral type need to be known. Since an homogeneous calibration 
extending from O stars down to early A stars is missing, we summarize all the adopted values. 



For O-type stars, we used the T* P ff from Marti ns et al.l (|2005l ). For early B supergiants, we 
adopted the T e g values given bv lCrowther et al.l (|2006bl ). while, for late B and A supergiants, those 



Table 7. List of field stars spectroscopically observed. 



ID 




Ra 






Dec 




B 


V 


R 


J 


H 


K s 


[3.6] 


[4.5] 


[5. 


8] 


[8.0] 


tel. 


Sp 


10 


18 


13 


47 


56 


-17 


57 


01 


43 


15 32 




13 64 


10 76 


9 


89 


9 60 


9 31 


9 34 


g 


22 


9 20 


ukirt 


K0 5III 


27 


18 


13 


14 


39 


-17 


54 


47 


87 






17 48 


9 50 


7 


84 


7 08 


6 90 


6 81 


Q 


55 


6 54 


uftirt- keck 




28 




13 


35 




-17 


54 


26 




12 85 


12 70 






7 


53 


7 21 


7 16 


7 22 


■j 


00 


7 03 




M4III 


29 


18 


13 


29 


14 


-17 


52 


27 


18 


16 19 


13 16 


13 13 


9 00 


7 


94 


7 50 


7 26 


7 51 


7 


26 


7 28 


keck 


M3III 


30 


18 


13 


26 


87 


-17 


55 


27 


50 








10.35 


8 


40 


7.52 


7.07 


7.12 


g 


90 


6.87 


keck 


M5III 


31 


18 


13 


16 


80 


-17 


55 


40 


02 


21 42 


16 78 


14 36 


9 79 


8 


47 


7 93 


7 63 


7 72 


7 


58 


7 55 


keck 


A/T9 "^TTT 


32 


18 


13 


16 


35 


-17 


50 


29 


89 








13 53 


9 


85 


8 08 


6 82 


6 80 


Q 


41 


6 47 


keck 


M2 51 


33 


18 


13 


25 


93 


-17 


57 


39 


84 


16 25 


16 37 




1161 


9 


27 


8 24 


7 55 


7 66 


7 


39 
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A/T 1 ^ "^TTT 


34 


18 


13 


28 
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50 
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26 
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57 


33 


25 
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9 
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9 


19 
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70 


-17 


51 


14 
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9 


76 


9.40 


9.21 


9.25 


9 


19 


9.24 


ukirt 
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57 


18 


13 


17 


23 


-17 


55 


21 


81 


15.93 


13.67 


15.02 


10.61 


9 


83 


9.56 


9.36 


9.45 


9 


33 


9.23 


ukirt 


K0.5III 


58 


18 


13 


12 


99 


-17 


54 


02 


32 








14.55 


11 


42 


10.02 


9.06 


9.01 


8 


71 


8.75 


ukirt 


> G 


59 


18 


13 


20 


01 


-17 


53 


50 


13 


15.83 


14.87 


15.93 


11.54 


10 


73 
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10.32 


10.42 


10 


34 


10.15 


ukirt 


K1.5III 


60 


18 


13 


20 


65 


-17 




50 
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12 


65 
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70 


9.70 
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31 
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00 
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Note. — For each star, number designations and coordinates (J2000) arc followed by magnitudes measured in different bands. J,H, and K s measure- 
ments arc from 2MASS, while the magnitudes at 3.6 /jm, 4.5 /^m, 5.8 /^m, and 8 (im are from GLIMPSE. B, V, and R associations arc taken from the 
astromctric catalog NOMAD. 
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1969 1. 



by IHumphrevs & McElrovl (| 19841 ). For B giants, we adopted the T P f? by Humphreys fe McElrov 



1984). For B and A dwarfs, we used T e g estimated by IHumphrevs & McElroyj (|1984r ) and 



Johnson 



For O-type stars, we used the bolometric corrections in K-band ( BCk) by iMartins fe Plez 



(|2006l ) . For early B-type stars, we used those provided by iBibby et al.1 (120081 ) . We estimated the 
BCk s of late B supergiants a nd gi ants by assuming bolometric cor rections in the V -ban d (BCy) 
from IHumphrevs fe McElrovl (|1984l ). and intrinsic V — K colors from iKoornneej (|1983l ) and IWegner 
( 19941 ). For dwarf stars, a set of homogeneous BCk was obt ained by interpolating an isochrone of 



0.5 Myr and solar metallicity f r omlLeieune fe Schaererl (|200ll ) at the assumed effective temperatures 
([Humphreys fe McElrovl 11984 : |johnsonlll966l ). 
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Table 8. Infrared colors and -BCk of supergiant stars. 



Sp. T eff BC K J-K H-K 



Reference 



03 


42551 
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-0.21 


-0.10 


04 


40702 
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-0.10 


05 


38520 


-4.40 


-0.21 
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06 
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-0.21 


-0.10 
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-0.10 
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Bl 


21500 
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-0.03 


B2 


18500 
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B3 
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B5 


13700 
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0.00 
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B8 
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0.12 
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Table 9. Infrared colors and BCk of giant stars. 



Sp. T 



eff 



BC K J-K H-K 



Reference 



03 


42942 


-4.85 


-0.21 


-0.10 


04 


41486 


-4.70 


-0.21 


-0.10 


05 


39507 


-4.54 


-0.21 


-0.10 


06 


36673 


-4.37 


-0.21 


-0.10 


07 


34638 


-4.19 


-0.21 


-0.10 


08 


32573 


-4.00 


-0.21 


-0.10 


09 


30737 


-3.80 


-0.21 


-0.10 


BO 


30300 


-3.59 


-0.18 


-0.09 


Bl 


21100 


-2.65 


-0.19 


-0.10 


B2 


18000 


-2.16 


-0.18 


-0.10 


B3 


17100 


-1.96 


-0.12 


-0.06 


B5 


16300 


-1.81 


-0.10 


-0.06 


B8 


12550 


-0.94 


-0.03 


-0.02 


B9 


11400 


-0.55 


-0.02 


-0.01 



Mar tins et al. 

Mar tins et al. 

Mar tins et al. 

Mar tins et al. 

Mar tins et al. 

Martins et al. 

Mar tins et al. 
Humphreys &z Mc Elrov 
Humphreys &: Mc Elrov 
Humphreys &: Mc Elrov 
Humphreys Mc Elrov 
Humphreys &z Mc Elrov 
Humphreys fc McElrov 
Humphreys & McElrov 



2005 ) .Martins fc Plez ('2006) 
2005 ) .Martins fc Plez (2006) 
2005 ) .Martins fc Plez (2006) 
2005 ) .Martins fc Plez (2006) 
2005 ) .Martins fc Plez (2006) 
2005 ) .Martins fc Plez (2006) 
2005 ) .Martins fc Plez (2006) 
1984 UWegner 1994) 
1984 UWegner 1994) 
1984 UWegner 1994) 
1984 ) . fWegner 1994) 
1984UWegner 1994) 
1984 ). (Wegner 1994 ) 
1984). (Wegner 1994 ) 
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Table 10. Infrared colors and BCk of dwarf stars. 



Sp. T eff BC K J-K H-K 



Reference 



03 


44616 


-4.76 


-0.19 


-0.04 


04 


43419 


-4.70 


-0.19 


-0.04 


05 


41540 


-4.62 


-0.18 


-0.04 


06 


38151 


-4.39 


-0.17 


-0.04 


07 


35531 


-4.15 


-0.17 


-0.04 


08 


33383 


-3.98 


-0.16 


-0.04 


09 


31524 


-3.86 


-0.15 


-0.04 


BO 


29600 


-3.72 


-0.15 


-0.04 


Bl 


24150 


-3.21 


-0.13 


-0.03 


B2 


19700 


-2.60 


-0.12 


-0.03 


B3 


18700 


-2.46 


-0.11 


-0.02 


B5 


13800 


-1.43 


-0.06 


-0.01 


B8 


12200 


-1.00 


-0.03 


-0.00 


B9 


10600 


-0.52 


-0.00 


0.00 


AO 


9850 


-0.25 


0.02 


0.01 


A2 


9120 


-0.03 


0.04 


0.02 


A5 


8260 


0.28 


0.08 


0.03 



Martins et al. C2005 
Martins et al. (2005 
Martins et al. (2005 
Martins et al. (2005 
Martins et al. (2005 
Martins et al. (2005 
Martins et al. (2005 
Humphreys k McElrov (1984 
Humphreys k McElrov (1984 
Humphreys k McElrov (1984 
Humphreys k McElrov (1984 
Johnson (1966 
Johnson (1966 
Johnson (1966 
Johnson (1966 
Johnson (1966 
Johnson (1966 



Leieune k Schaerer 



Leieune k Schaerer 



Leieune k Schaerer 



Leieune k Schaerer 



Leieune k Schaerer 



Leieune k Schaerer 



Leieune k Schaerer 



Leieune k Schaerer 



Leieune k Schaerer 



Leieune k Schaerer 



Leieune k Schaerer 



Leieune k Schaerer 



Leieune k Schaerer 



Leieune k Schaerer 



Leieune k Schaerer 



Leieune k Schaerer 



Leieune k Schaerer 



(2001 ) 
(2001) 
(2001) 
(2001) 
(2001) 
(2001 ) 
(2001 ) 

(2001) 
(2001) 
(2001) 
(2001) 
(2001) 
(2001 ) 
(2001 ) 
(2001 ) 
(2001 ) 
(2001) 



Table 11. Infrared colors and BCk adopted for RSGs and WRs. 



Sp. 


T eff 


BC^ 


J-K 


H-K 




Reference 


K2I 


4015 


+2.5 


+0.65 


+0.13 


LevesQue et al. 


(2005^1 .Koornneef (1983) 


Mil 


3745 


+2.7 


+1.00 


+0.22 


Levesaue et al. 


(2005).Koornneef (1983) 


M2.5I 


3615 


+2.8 


+1.06 


+0.27 


Levesaue et al. 


(20051 .Koornneef (1983) 


M3.5I 


3550 


+2.9 


+1.16 


+0.28 


Levesaue et al. 


(20051 .Koornneef (1983) 


WN7o 




-3.9 


+0.13 


+0.11 




Crowther et al. (2006a) 


WN7b 




-3.5 


+0.37 


+0.27 




Crowther et al. (2006a) 



